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ABSTRACT 

The uncertain origin of the recently-discovered ‘changing-looking’ quasar phenomenon - in which a 
luminous quasar dims significantly to a quiescent state in repeat spectroscopy over ^10 year timescales 
- may present unexpected challenges to our understanding of quasar accretion. To better understand 
this phenomenon, we take a first step to building a sample of changing-look quasars with a systematic 
but simple archival search for these objects in the Sloan Digital Sky Survey Data Release 12. By lever¬ 
aging the >10 year baselines for objects with repeat spectroscopy, we uncover two new changing-look 
quasars, and a third discovered previously. Decomposition of the multi-epoch spectra and analysis 
of the broad emission lines suggest that the quasar accretion disk emission dims due to rapidly de¬ 
creasing accretion rates (by factors of >2.5), while disfavoring changes in intrinsic dust extinction 
for the two objects where these analyses are possible. Broad emission line energetics also support 
intrinsic dimming of quasar emission as the origin for this phenomenon rather than transient tidal 
disruption events or supernovae. Although our search criteria included quasars at all redshifts and 
transitions from either quasar-like to galaxy-like states or the reverse, all of the clear cases of changing- 
look quasars discovered were at relatively low-redshift {z ^ 0.2 — 0.3) and only exhibit quasar-like to 
galaxy-like transitions. 

Subject headings: galaxies: active, quasars: emission lines, quasars: general 


1. INTRODUCTION 

The quasar phenomenon is thought to be a relatively 
brief stage of galaxy evolution involving rapid accretion 
onto the central supermassive black hole (S MBH; e.g., 
Salpeter||I964| |Lynden-Bell||1969t |Rees||I978|) . Observa- 
tional constraints on litetirnes show that quasar p hases 
in galaxies generally last for a total of I0^~ ^ years ( Mar- 
tini & Weinberg 2001 Kelly et al. 20101, after which 
the accretion rate drops dramatically and the active nu¬ 
cleus transitions to a low-luminosity ac tive galactic nu¬ 
cleus (AGN) or quiescent galaxy state (|Churazov et al. 
20051. Cosmological simulations of galaxy formation 
that include sub-grid models of SMBH growth and feed- 
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back have suggested that the accretion history of SMBHs 
may be episodic, where luminous quasar p hases are reg¬ 
ulated by quasar feedback processes (e.g.,|Di Matt eo et 
al. 2005 iHopkins et al. 2005 Springel et aT 2005). Al¬ 
thoughtire exact characteristics of quasar light curves 
over cosmic time are difficult to infer observa tionally, in¬ 
direct argu ments based on AGN populations (iSch awinsMi 
2015 1 and properties of their host-galaxies (|HicEox| 


et al. 


et al. 


20141 have also suggested that AGNs dramatically 


‘dicker’ in luminosity between luminous quasar and qui¬ 
escent galaxy phases. However, direct observations of 
such transitions in luminous quasars have thus far been 
scarce. 

It has been suggested that the transition from quasars 
to low-luminosity AGN or quiescent galaxies may not 
be observable in individual objects due to the long 
timescales expected for this process. Such dramatic 
changes in the accretion state are commonly observed 
in X-ray binaries, which can undergo spectral state 
transitions between the high-luminosity/soft-spectrum 
and l ow-luminosity /hard-spect rum states in the X-rays 
(e.g., Homan & Belloni 2005). Scaling the hours-long 
timescales observed tor these spectral state transitions in 
X-ray binaries to ^10® Mq SMBHs predicts transition 
times cales in quasars of ~10*“® years (Sobolewska et al. 


2011). Indirect evidence for a luminous quasar transi- 


tion m an individual object was previously provided by 
observations of Hanny’s Voorwerp, a serendipitously dis¬ 
covered ionized emission-line gas cloud l ying ~20 kpc 


away from the quiescent galaxy IC2497 (Lintott et al. 


2009). Based on multi-wavelength observations, it is ar- 
gued that this gas cloud could only have been ionized by 
an AGN continuum with luminosity >10"*® erg s“*; this 
implies that the nearby quiescent galaxy was recently in 
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a quasar state, with a transition timescale of years 
and consistent with expectations (Schawinski et al.|2010 


Keel et al.|2012a). Since this discovery, many other can- 


didate lading AGNs with extended emission-line regions 
have been found and investigate d, resulting in similar in - 
ferred transition timescales (e.g.. Keel et al. 2012b 2015). 

Previously, transitions of Seyfert 1 galaxies to Seyfert 
1.8/1.9 (and/or vice versa) have been directly observed in 
repeat optical spectroscopy over timescales of ^10 year s 


>P ^ 

dGoodrich l lOOSi fShappee et al. 2014 Denney et al. 2014|). 

Goodrich| (|1995p showed that the origin of the observed 
transitiori in some of these AGNs are consistent with in¬ 
trinsic changes in the AGN continuum emission, while 
variations in dust obscuration along the line of sight is 
favored for others, although these two effects may occur 
in co ncert if dust is embed ded in the narrow-line region 
gas (Netzer & Laor 1993). Similar behavior has also 
been previously observed in the X-ray spectra of Seyfert 
galaxies (termed “changing-look AGN”), also interpreted 
as being due to either dramatic changes in obs curation or 
intrinsic changes in the nuclear emission (e.g. 

20021 Matt et al.||2003[ IPuccetti et al.||200R 


Guainazzi 
Risaliti et 


^ 2009 Marchese et al.||2012p. The recent discovery of 


the hr st changmg-look (GL) :(uasar by |LaMassa et al. 
(2015) extends this transitional phenonienon to AGiNs 
in new luminosity a nd redshift regimes (see Figure 1 of 
LaMassa et al.|2015 ). Repeat optical spectroscopy of this 
luminous quasar (SDSS J015957.64-|-(103310.5, hereafter 
J0159-I-0033) shows a dramatic decrease in the quasar 
continuum emission, accompanied by disappearance of 
the broad H/3 line and strong dimming of the broad 
Ha line, showing that a simple orientation-based view 
of AGN unification is incomplete. Surprisingly, the ob¬ 
served transition in this GL quasar occurred over rest- 
frame timescales of ~7 years; this is much shorter than 
the ^10"* year timescales expected for this transition to 
occur from previous arguments based on X-ray binaries 
and extended emission line regions surrounding quiescent 
galaxies. 

The origin of chan ging-look behavior in luminous 


quasars is uncertain. LaMassa et al. (2015) demon¬ 


strated that the observed dimming of the quasar con¬ 
tinuum in J0159-I-0033 coincides with broadening of the 
broad Balmer emission lines, such that the derived black 
hole mass (estimated through single-epoch spectroscopic 
black hole mass methods) is preserved. This behavior 
is consistent with intrinsic dimming of the quasar con¬ 
tinuum emission, while a scenario in which the contin¬ 
uum and emission line dimming is caused by an increase 
in dust extinction is disfavored t hrough modeling of the 
spectral changes. Furthermore, LaMassa et al. (2015) 
also argue that obscuration by a dust cloud outside the 
broad line region in a circular Keplerian orbit is unlikely 
since its crossing time across the broad line region would 
be much longer than the observed transition timescale. 

If the dimming of quasar emission in CL quasars is 
intrinsic, then the observed behavior can be caused by 
dramatic changes in the accretion flow, which may oc¬ 
cur during transitions betwe en radiatively efficient and 
inefficient accretion regime s (|lchimaru 1977 Rees et al. 
1982 Narayan & Yi 1994). Ihermal and dynamical in¬ 
stabilities in the accretion disk may also produce strong 
changes in the disk emission on even shorter tim e scales 
(Lin & Shields 1986 Siemiginowska et al. 1996). Mer¬ 


loni et al. (2015) argued that the LaMassa et al. (2015) 
GL quasar may instead be a transient stellar tidal dis¬ 
ruption event (TDE) near the central SMBH (see also 
Eracleous et al. 1995), which would cause a luminous 
nuclear flare, followed by a slow dimming over the few- 
years timescales observed. This scenario is supported by 
their image-differencing light curves of this GL quasar, 
which show that the time evolution of the broadband 
nuclear emission is consistent with that expected from 
TDEs. In any of these scenarios, GL quasars represent 
an intriguing new phenomenon that can provide unique 
insights to AGN accretion and structure, and warrant 
additional investigation. 


Since the discovery of LaMassa et al. (2015), at least 
two more spectro scopic Ol quasars h ave been serendipi- 
tously discovered (Runnoe et al. 2016), primarily through 
visual inspection ot early spectroscopic data from the 


Time -Domain Spectroscopic Survey (Morganson et al. 
2015) in the Sloan Digital Sky Survey IV (SDSS IV). 


These results suggest that CL quasars may be surpris¬ 
ingly common, and can be found by mining spectral data 
sets with large numbers of repeat quasar and galaxy 
spectra over a sufficiently long baseline. Motivated by 
these results, we present a simple archival search for 
CL quasars in SDSS-I/B/BI Data Release 12 (DR12, 
Alam et al. 2015), which includes a total of ~4 x 10® 
optical/near-IR spectra over ^10,000 deg^ of sky. This 
data set includes a significant number of repeat spectra 
spanning a >10 year baseline, and many of its various sci¬ 
ence programs specifically targeted quasars and galaxies. 
We aim to cast a wide net in this investigation, and in¬ 
clude quasar and galaxies at all redshifts while remaining 
impartial in our search with regard to whether each ob¬ 
ject transitions from a quasar-like to galaxy-like state or 
vice versa; this approach could not only yield a sample 
of CL quasars, but also potentially a spectroscopic tran¬ 
sition from quiescent galaxy to luminous quasar, which 
would have intriguing implications for their origin. 

The outline of this paper is as follows: Section 2 
presents the data sets and criteria used in our search. 
In Section 3, we describe the changing look quasars 
found in our search and modeling of their broad emis¬ 
sion lines. Section 4 evaluates evidence from our sample 
favoring various scenarios for the origin of CL quasars. 
We summarize and conclude in Section 5. Throughout 
this paper, we assume a standard ACDM cosmology with 
Dm = 0.309, Da = 0.691, and Hq = 67.7 km s"! Mpc~\ 
consistent with the Planck full-mission results of IPlanckI 
Collaboration et al. (2015). 


2. AN ARCHIVAL SPECTROSCOPIC SEARCH 
2.1. Search Criteria 

We utilize the list of all 4,355,202 spectra in SDSS 
( York et ar]|2000 ) DR12, and perform the selection cuts 
detailed below to produce a final sample of 117 CL quasar 
candidates. These spectra were taken by the SDSS 2.5m 
telescope (Gunn et al. 2006) using the SDSS-I/ H and 
Baryon Oscillation Spectroscopic Survey (BOSS, Eisen-| 


stein et al.||201ll 

([Smee et al.||2Ui; 


Dawson et al. 


2013) spectrographs 
2013p, and compiled m the ‘spAlT files 
produced by the SP SS spectroscopic reduction pipeline 
(Bolton et al. 2012). Eor a CL quasar to show a con¬ 


vincing transition, its multi-epoch spectra must clearly 
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Table 1 

Measured SDSS spectral properties of the changing-look quasars in our sample. 


SDSS 

Object 

z 

MJD 

Ho FWHM’' 
[km s“*-] 

logioLua'' 

[erg 8“*] 

H^ FWHM® 
[km s“*] 

logioLuy 

[erg s”*] 

logioALsioo 
[erg s”*] 

J015957.64-I-003310.5*’ 

0.312 

51871 

3788 ± 163 

42.36 ± 0.04 

4714 ± 682 

41.88 ± 0.63 

43.52 ± 0.05 



55201 

5954 ± 857 

41.72 ± 0.11 


<41.32® 

43.27 ± 0.06 

J012648.08-083948.0 

0.198 

52163 

4121 ± 223 

42.00 ± 0.04 

4297 ± 1165 

41.55 ± 0.20 

43.43 ± 0.03 



54465 


<40.20'*, <40.51® 


<40.17*", <40.28® 

<42.30® 

J233602.98-1-001728.7 

0.243 

52096*' 

6289 ± 1180 

41.86 ± 0.20 

6993 ± 2271 

41.28 ± 0.20 

43.04 ± 0.09 



55449 

7209 ± 1367 

41.48 ± 0.22 


<40.60® 

42.56 ± 0.18 


^ These measurements of the luminosities a nd widths are for the b road components of these Balmer lines. 

^ Changing-look quasar previously found by|LaMassa et al.|(|2015ll and also disused in|Merloni et aL|||2015|l. 

^ 5(T upper limit assuming FWHM predicted trom Ma, see Section 3.2. 

^ 5(7 upper limit assuming FWMH of 7960 km s~^, see Section 3.2. 

® 5(7 upper limit assuming FWMH of 4121 km s~^, see Section 3.2. 

^ 5(7 upper limit assuming FWMH of 8920 km s~^, see Section 3.2. 
s 5(7 upper limit assuming FWMH of 4297 km s“^, see Section 3.2. 

^ This MJD is the mean of four closely-spaced epochs of spectra that have been stacked (see discussion in Section 3.1). 


Table 2 

Inferred SDSS spectral properties of the changing-look quasars in our sample. 


SDSS 

Object 

MJD 

logloAYsH.Ha 

logioMBH,H/9 

Wq] 

logic (Lbol/ABdd,Ha) 

logic (Abol / Afidd,H/3 ) 

J015957.64-I-003310.5® 

51871 

55201 

7.93 ± 0.10 

8.20 ± 0.26 

8.02 ± 0.33 

-1.6±0.1 
-2.1 ± 1.6 

-1.7± 1.3 

J012648.08-083948.0 

52163 

54465 

7.96 ± 0.10 

7.89 ± 0.84 

-1.7± 0.4 

-1.7± 0.7 

J233602.98+001728.7 

52096*' 

55449 

8.13 ± 0.29 

8.00 ± 0.30 

8.11 ± 0.53 

-2.3 ±0.4 
-2.7 ± 2.6 

-2.3 ±0.7 


^ Changing-look quasar previously found by |LaMassa et ar] l |201^ and discussed in [Merloni et al.| ( [2015^ . 

This MJD is the mean of four closely-spaced epochs of spectra that have been stacked (see discussion in Section 3.1). 


possess quasar-like spectral features in one epoch (power- 
law continuum and broad emission lines), and galaxy-like 
features in another epoch (absorption spectra and nar¬ 
row emission lines if star-formation or nuclear activity 
is present). While a sophisticated method of detecting 
this transition in repeat spectra of each object is likely 
to be more sensitive to subtle changes, our current goal 
is to search only for the most obvious and convincing 
cases of CL quasars. Thus, our simple approach relies 
on the automated SDSS pipeline to classify each spec¬ 
trum as quasar-like or galaxy-like. Specihcally, using 
the CLASS spectral classification provided for each spec¬ 
trum in our sample (which is based on fitting to a set 


of galaxy , quasar, and stellar eigenspectra, see [Bolton et 
al. 20121, we create two subsamples: a galaxy-like sam¬ 
ple of 2,510,060 spectra where CLASS = ‘GALAXY’, and 
a quasar-like sample of 587,306 spectra where CLASS = 
‘ QSO ’. In both these subsamples, sky hbers have been 
removed using the sourcetype targeting keyword. Al¬ 
though it is well-known that these automated pipeline 
classifications occasionally fail to accurately classify the 
observed spectra, we emphasize that we are focusing 
on the most obvious and convincing CL quasars, for 
which the pipeline classifications will suffice for this ini¬ 
tial archival search. 

For each spectrum in our galaxy-like sample, we po¬ 


sitionally match to the quasar-like sample using a 1" 
matching radius to identify objects with repeat spec¬ 
tra and disparate classifications. We impose additional 
quality-control conditions on this search: the difference 
in the pipeline redshifts between the two epochs must 
be |Az| < 0.01, and the absolute value of the rest-frame 
time lag between the two epochs is |Afrest-framel > 4 
years. These conditions remove a signihcant number of 
false-positives in which the pipeline catastrophically fails 
to fit the spectrum in one of the epochs. This search re¬ 
sults in 180 pairs of repeat spectra of 117 unique objects 
(a few objects have more than two epochs of spectra), 
and we visually inspect all spectra of each of these CL 
quasar candidates. 

From the visual inspection, we find three clear cases 
of CL quasars ( listed in Table 1), which include the 


LaMassa et al. (2015) CL quasar (J01595-I-0033) and 
two additional new convincing cases (J01264—0839 and 
J2336-I-00172). All of these CL quasars exhibit quasar¬ 
like to galaxy-like transitions, and are at relatively low 
redshifts of z ~ 0.2 — 0.3 in comparison to the parent 
sample of SDSS quasars, wh ich are overwhelm ingly at 
z > 1 (e.g. see Figure 2 in Paris et al. 2014). In all 
three of the CL quasars discovered in our search, the 
broad H/3 emission disappears while the broad Ha emis¬ 
sion dims significantly (and disappears in J01264—0839), 
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accompanied by dimming of blue quasar UV/optical con¬ 
tinuum emission. The median seeing during the expo¬ 
sures of these three CL quasar candidates ranged from 
1.43"to 1.93", consistent with typical values for SDSS 
spectra (i.e. these exposures are not strong outliers). 
Small differences in the seeing between exposures at this 
level are not responsible for the observed spectral vari¬ 
ability (including broadening of the broad emission lines, 
see Section 3.2). We additionally find one ambiguous CL 
quasar where the latest epoch of SDSS spectra appears 
to show a galaxy-like spectrum at blue wavelengths, but 
is corrupted at the redder wavelengths. In Appendix 
A, we present additional recent non-SDSS spectroscopy 
which demonstrates that this object does not transition 
to a galaxy-like state at the epoch of the latest spectrum; 
this behavior instead likely stems from known instrumen¬ 
tal issues affecting the particular fiber of the corrupted 
SDSS spectrum. 

In the visual inspection of the multi-epoch spectra, 
the vast majority of the false-positives from our search 
were cases where the pipeline switched between CLASS = 
‘GALAXY’ and CLASS = ‘QSO’ classifications in repeat 
spectra despite little change in the spectral properties. 
Often, this occurs for AGN at redshifts of z ~ 0.4 in 
which the broad Ha emission line is redward of the 
smaller wavelength coverage of the SDSS spectrograph in 
the earlier epoch (leading to a CLASS = ‘GALAXY’ clas¬ 
sification), but visible in the later epoch from the BOSS 
spectrograph due to its slightly larger wavelength cover¬ 
age (leading to a CLASS = ‘QSD’ classification). As part 
of the visual inspection, the fiber plugging positions and 
targeting flags for each pair of repeat spectra were com¬ 
pared to ensure that there is no offset in the fiber posi¬ 
tion between the two spectra, artificially leading to more 
host-galaxy emission in the SDSS spectrum. This offset 
can occur since some SDSS fibers were part of a SDSS 
program to test the redshift recovery of the spectroscopic 
pipeline in SDSS-III relative to that in SDSS-I/II, and 
are identified using the PROGRAM = ‘APB I AS’ target flag 
in the spectra as well as their disparate fiber plugging 
positions in repeat spectra. One of these objects was re¬ 
covered in our search, and was removed from our sample. 
Aside from this offset object and the three CL quasars, all 
the remaining 113 candidates were rejected because our 
visual inspection of their repeat spectra did not reveal 
dramatic disappearance or appearance of broad emission 
lines. In Appendix B, we further discuss and show exam¬ 
ples of the CL quasar candidates rejected in our visual 
inspection. 

3. SPECTRAL PROPERTIES OF 
CHANGING-LOOK QUASARS 

In this section, we describe our decomposition of the 
multi-epoch spectra of our three CL quasars into host- 
galaxy and quasar components. We then fit the broad 
emission lines in the quasar components and analyze 
their properties. Although single-epoch broadband imag¬ 
ing of our CL quasars is available through SDSS, we do 
not compare the multi-epoch spectra to their photometry 
in our analysis. This is primarily because our CL quasars 
are spatially extended (i.e. resolved) in the imaging, and 
the SDSS fibers may have different diameters depend¬ 
ing on epoch. This causes the stellar contamination in 
the imaging photometry to vary with aperture, compli¬ 


cating comparisons between the imaging and the multi¬ 
epoch spectra. Our approach to investigating the nuclear 
emission of our CL quasars instead relies on decomposi¬ 
tion of the spectra to remov e contamination from the 
host galaxy. We note that the LaMassa et al. (2015) CL 
quasar recovered in our search lies in the SDSS Stripe 82 
footprint; image differencing light curves of its nuclear 
emissi on are presented and discussed in Merloni et al. 


3.1. Spectral Decomposition 

We decompose the quasar and host-galaxy components 
in the spectra from both epochs for each of the three 
CL quasars found in our search. For the spe ctral de¬ 
composition, we follo w the gener a l met hod of Vanden 
Berk et al. (2006) and Shen et al. (2015), with only mi¬ 
nor modifications, which is based on htting quasar and 
galaxy eigenspectra. Specifically, the CL quasar spec¬ 
tra are each fitted using a mix of quasar and galaxy 
eigenspectra that are created from a principal compo¬ 
nent analysis of large samples of SDSS quasar and galaxy 
spectra. This decompositio n method using e igensp ectra 
differs from the approach of [LaMassa et al. (2015), who 
instead fit a model of power-law quasar continuum emis¬ 
sion, host-galaxy emission generated from a stellar pop¬ 
ulation synthesis model, and emission lines to the multi¬ 
epoch spectra. The main advantage of our approach is 
that the continuum emission is fitted and decomposed 
empirically without having to rely on the accuracy of a 
power-law and stellar population synthesis models. For 
example, a mixture of eigenspectra can more accurately 
describe a galaxy spectrum consisting of stars with a con¬ 
tinuous range of ages in comparison to the simple star for¬ 
mation histories assumed by stellar population synthesis 
models. The primary disadvantage of our approach is 
that because narrow emission lines are present in both 
the quasar and galaxy eigenspectra, separating the nar¬ 
row line emission from the continuum emission is less 
straightforward (this does not adversely affect the broad 
emission lines). Although the CL quasar J0159-I-0033 
recovered in our search was previ ously discovered and 
analyzed by [LaMassa et al. (2015), we nevertheless in¬ 
clude it in our analysis below to demonstrate whether 
our independent spectral decomposition and broad emis¬ 
sion line fitting for this object produces results that are 
consistent with their published values. 

We first correct all spectra in our sample for Galac¬ 
tic ext inction, using the maps of |Schlafiy fc Finkbeiner 
(12011 ) an d the Milky Way reddening law of Cardelli et 
(|l989). To facilitate the spectral decomposition, we 
resample all our spectra and the eigenspectra to a com¬ 
mon wavelength grid of the form logioA = 3.35 -I- 0.001a, 
for integer a from 0 to 5,914. The wavelength cover¬ 
age of this common grid is wide enough to accommo¬ 
date all spectra in our sample, and is similar to the na¬ 
tive SDSS resolution. For our spectral decomposition, 
we utilize the eigenspectra from the principal c omponent 
analysi s of ~ 17,000 SDSS quasar spectra from | Yip et al. 
(|2004a|, a nd 170,000 SDSS galaxy spectra from | Yip et 
al.|(|2004b|), respectively. Specifically, we fit combinations 
of the hrst five quasar and first five galaxy eigenspectra, 
with their amplitudes (i.e., PCA coefficients) as 10 free 
parameters. The se fits are performe d through a simple 
minimization. Yip et al. ( 2004a|b ) demonstrated that 
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Figure 1. Top: Spectral decomposition for the two epochs of spectra of SDSS J015957.64-1-003310.5 (see|LaMassa et al.||2015|l, in the 
rest-frame. The black lines are the observed spectra, and the green and blue lines are the reconstructed quasar and host-galaxy spectra 
from the eigenspectra decomposition, respectively. The best-fit model to the observed spectrum from the decomposition (i.e. sum of the 
green and blue lines) is the red line. The dramatic dimming in the quasar continuum and broad Balmer emission lines are consistent with 
intrinsic dimming of the accretion disk emission rather than dust extinction. Bottom left: Fitting of the H/3 line region in the decomposed 
quasar spectrum, for the two epochs of spectra. The decomposed quasar spectra are the black lines, the best-fit broad and narrow H/5 
emission lines are the blue lines, and the total fits to the decomposed quasar spectra (including quasar continuum and all emission lines) 
are shown in red. Although narrow emission lines are included in the fit, their amplitudes in the decomposed spectrum are not equivalent 
to the narrow emission lines in the observed spectrum since they are partially subtracted as part of the host-galaxy spectrum (see Section 
3.2). Bottom right: Similar to the bottom left panels, but for the Ha emission lines. 
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Figure 2. Spectral decomposition for the two epochs of spectra of SDSS J012648.08-083948.0 in the rest-frame, same format as Figure 1. 


the first five eigenspectra in their PCA analysis captured 
98.29% and 98.37% of the variance in their quasar and 
galaxy spectra samples, respectively. Not surprisingly, 
we hnd that extending the spectral decomposition to the 
first 10 quasar and galaxy eigenspectra did not notice¬ 
ably improve the resulting fits to the observed spectra, 
and so we only utilize the first five in our analysis for 
simplicity. 

Since the Petrosian radius measured in SDSS imaging 
for our three CL quasars are between 1.45" to 2.58", they 


are spatially extended (i.e. resolved). The host-galaxy 
and quasar contributions in each resulting SDSS spec¬ 
trum are thus dependent on the fiber diameter. SDSS- 
I/II spectra were obtained using 3" diameter fibers, while 
SDSS-III (BOSS) spectra were acquired using 2" diam¬ 
eter fibers. For J0159-I-0033, we decompose the spec¬ 
tra from the two epochs separately because the earlier 
spectrum was obtained with a 3" fiber, while the later 
spectrum was obtained with a 2" hber. In contrast, 
both spectra of J0126—0839 were obtained with 3" hbers, 
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SDSSJ233602.98+001728.7, z = 0.243 



rest-frame wavelength [A] 

H/3 Ha 



Figure 3. Spectral decomposition for the two epochs of spectra of SDSS J233602.98-1-001728.7 in the rest-frame, same format as Figure 
1. We note that the earlier epoch of spectra at MJD 52096 is a mean stack of four epochs of SDSS spectra taken within a 2-year period 
(MJD of 51783, 51877, 52199, 52525), during which no strong spectral changes were observed. 


therefore the host-galaxy contribution should be constant 
between the two spectra. For this object, we decom¬ 
pose both epochs of spectra simultaneously and impose 
the additional constraint of constant galaxy parameters 
between the spectra from the two epochs. Finally, for 
J2336+0017, a total of five epochs of SDSS spectra are 
available, including four early epochs (with 3" fibers) 
within a 2-year timespan in the observed frame during its 


quasar-like phase (MJD of 51783, 51877, 52199, 52525), 
and an epoch 9 years later with a 2" fiber. Since the 
spectral changes in the four early epochs of this object 
are relatively small, we simply use the mean spectrum 
of these four epochs in our spectral decomposition to 
achieve higher signal-to-noise, and our quoted MJD of 
52096 is thus actually the mean MJD of the four early 
epochs. To accommodate the different fiber diameter of 
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the early and later spectra, the galaxy parameters in the 
fitting are allowed to vary between the mean earlier spec¬ 
trum and the later spectrum. In Section 4.2, we discuss 
the evolution in the continuum luminosity of this object 
over its five separate spectroscopic epochs. 

The results of our spectral decomposition for each of 
the three CL quasars are shown in Figures 1-3, which dis¬ 
play both observed epochs of each CL quasar, as well as 
their best-fit quasar and galaxy components from eigen- 
spectra. Following previous conventions, we refer to the 
fitted quasar and galaxy spectra shown in Figures 1-3 as 
the ‘reconstructed’ quasar and galaxy spectra, while the 
‘decomposed’ quasar spectra (not shown) are the original 
spectra with their corresponding reconstructed galaxy 
spectra subtracted. All of our broad emission line analy¬ 
sis in Section 4.2 is performed on the decomposed quasar 
spectra rather than the reconstructed quasar spectra, 
since this allows us to use the uncertainties on the flux 
densities from the original spectrum. 


3.2. Broad Emission Line Analysis 

Using our decomposed quasar spectra in both epochs 
for each of our CL quasars, we measure the proper¬ 
ties of the Ha and H/1 broad emission lines to estimate 
black hole masses Mbh and bolometric Eddington ratios 
Lboi/AEdd- Our broad emission line fitting proced ure 
generally follows the method of Shen et al. (20111, in 
which the single-epoch virial Mbh estimates are based 
on the broad emission line Full-Width Half Maximum 
(FWHMs), as well as a radius-luminosity relation for 
the broad line region from reverberation mapping of low- 
redshift AGNs. We emphasize that the narrow emission 
lines in the decomposed quasar spectra (observed spec¬ 
trum minus the best-fit galaxy spectrum) we use for the 
broad emission line fitting are not equivalent to those in 
the original observed spectrum. This is because narrow 
emission lines are also present in the galaxy eigenspec- 
tra (and thus the best-ht host-galaxy spectrum), which 
is subtracted to obtain the decomposed quasar spectra. 
However, we include the narrow emission lines in our 
analysis below to avoid biases in fitting the broad emis¬ 
sion lines. 

For the Ha region of each object, we use the decom¬ 
posed quasar spectrum and fit the local continuum emis¬ 
sion in the wavelength windows of [6400, 6500]A and 
[6800, 7000]A to a power-law. In the Ha line wavelength 
window of [6500, 6800]A, we fit for the narrow Ha com¬ 
ponent, the [N II] AA6548,6584 doublet, and the [S H] 
AA6717,6731 doublet, using a single Gaussian for each 
emission line. The redshifts of the narrow lines are con¬ 
strained to be the same, and their widths are constrained 
to be <1200 km s“^. The broad component of the Ha 
emission is fit with a Gaussian with width constrained to 
be >1200 km s“^, with its central wavelength as a free 
parameter. 

Similarly, for the H/3 region, we fit a local power-law to 
the continuum wavelength windows of [4435, 4700]A and 
[5100, 5535]A, and we fit emission lines in the H/3 line 
wavelength window of [4700, 5007]A. In the continuum 
w avelength window, we incl ude the optical Fe H template 
of Boroson & Green (1992) in the fit. However, for each 
of our three GL quasars m the dimmer galaxy-like epoch 
and two in the brighter quasar-like epoch (J012648.08- 


083948.0 and J233602.98-1-001728.7) , the Fe H emission 
is weak and the template fit is poorly constrained. Thus, 
the Fe H template is included in the continuum fit only 
for J0159-I-0033 in its brighter quasar-like spectral epoch. 
For the narrow lines, we fit single Gaussians for the nar¬ 
row component of H/3 and the [O HI] AA4959,5007 dou¬ 
blet, with widths constrained to be <1200 km s“^ and 
redshifts constrained to be the same. The broad H/3 com¬ 
ponent is fit with a single Gaussian with width >1200 km 
s“^, and its central wavelength is a free parameter. 

The above spectral fitting produced the FWHM and 
luminosities of the broad Ha and H/3 components, as well 
as the quasar continuum luminosity at 5100A AL 5100 , for 
each epoch of our GL quasars. These measured proper¬ 
ties of the broad emission lines are tabulated in Table 
1, and the spectral fits are also presented in Figures 1- 
3. All uncertainties are calculated through 10^ Monte 
Garlo realizations of each spectrum based on their Icr flux 
density uncertainties. The spectral htting procedure is 
performed for each resampled spectrum, and the uncer¬ 
tainties quoted for each parameter are thus the Icr spread 
in the resulting distributions of resampled parameters. 

Using the Ha broad emission line FWHM and lumi¬ 
nosities from our fits (where detected), single-epoch black 
hole masses Mbh Hr> are estimated using the relation 
from Greene et al. ( |2010 ) of 


MBH.Hct = 9.7 X 10® 


AL, 


FWHM(Ha) 
1000 km s“i 


-I 2.06 


1 0.519 


5100 


10^^^ erg s 


-1 


M 
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( 1 ) 


based on the radius-luminosity relation of |Bentz et al.| 
(2009). To verify these Ha-based estimates, we also es- 


timate single-epoch black hole masses based on our H/3 
broad emission line fits Mbh.h /i using the relation from 
Vestergaard & Peterson (2006) of 


.n6.9i rFWHM(H/3) 
[lOOOkms-i 
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^^5100 


10^^^ erg s 
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0.5 


M 


0- 
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For J0159-I-0033 and J2336-I-0017, the broad H/3 com¬ 
ponent is not detected in the fainter galaxy-like spec¬ 
trum. We provide 3cr upper limits on the broad H/3 lu¬ 
minosities for these two objects in Table 1. These upper 
limits are estimated in the fainter galaxy-like epoch by 
refitting the Gaussian broad H/3 line in 10® Monte Garlo 
resamplings of the spectrum, with the FWHM of the 
broad H/3 component fixed to the width predicted from 
the observable Ha FWHM using Equation 9 from |Shen| 

eTaTl (I20TTI). 

Tot JU12o—0839, neither broad Balmer emission nor 
quasar continuum emission are detected in the fainter 
galaxy-like spectral epoch, and thus the origin of the dim¬ 
ming in this object is poorly-constrained. We provide a 
5(7 upper limit on the continuum luminosity AL 5100 in 
Table I, through 10® Monte Garlo resamplings of the 
spectrum. To obtain limits for the luminosities of the 
Balmer emission lines, we must first assume a FWHM. 
Therefore, we report here the results for two different as¬ 
sumptions: (a) the FWHM in the dim state is the same 
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1233602.98 + 001728.7, Ha 



Figure 4. Decomposed quasar spectrum of J2336H-0017 in the Ha 
region from the earlier, quasar-like epoch (red), compared to the 
dereddened decomposed quasar spectrum from the later, galaxy¬ 
like epoch (blue). Although changes in dust extinction can rea¬ 
sonably explain the dimming of the continuum emission in this 
changing-look quasar, the extinction required cannot explain the 
strong changes in the broad emission line (see discussion in Section 
4.1), disfavoring an extinction origin for changing-look quasars. 


as that in the bright state, and (b) the FWHM increases 
as the continuum luminosity drops in order to preserve 
the black hole mass, as prescribed by Equations 1 and 2 
above. Since the continuum is not detected in the dim 
state, we can only obtain upper limits on the FWHM of 
Ha and H/3 of 7960 and 8920 km s“^, respectively. Thus, 
we assume these (fixed) values of the FWHM in order to 
obtain upper limits on the line luminosities. The limits 
we obtain for cases (a) and (b) are given in Table 1. We 
emphasize that since neither the continuum nor the lines 
are detected in the dim state, we cannot constrain the 
dimming mechanism for this object as well as did for the 
other two objects. Therefore, neither the intrinsic dim¬ 
ming nor the reddening hypothesis can be excluded for 
J0126-0839. 


We estimate the bolometric luminosity Lboi of our 
quasars by multiplying the continuum luminosity mea¬ 
sured from the decomposed quasar spectra at 5100A, 
AT,t;inn, by the bolom etric correction factor of 8.1 from 


Runnoe et al. (20121. The Eddington ratio is then 
Aboi/T^Edd = iboi/(l-3 X 10 ^®Mbh), for Mbh in units of 
Mq, and Lboi in units of erg s“^. These inferred quan¬ 
tities are tabulated for both epochs of our CL quasars 
in Table 2. The Mbh and Lboi/TEdd values derived for 
each epoch of spectra from Ha and H/? are consistent to 
within the uncertainties, validating our broad emission 
line fits. 


4. DISCUSSION 

In this section, we investigate the origin of the CL 
quasar phenomenon using our sample. In particular, we 
will focus on evidence for and against variable dust ob¬ 
scuration, tidal disruption events, Type Hn supernovae, 
and intrinsic dimming of quasar emission, based on a 
variety of approaches. We conclude from these investi¬ 
gations that CL quasars are likely to be due to intrinsic 
dimming of the nuclear emission, as a result of rapidly 


decreasing accretion rates. 


4.1. A Dust Extinction Origin? 

We assess whether extinction by an intervening dust 
cloud along the line of sight outside the broad line re¬ 
gion can cause the dimming in the continuum and broad 
emission li nes observed in our CL quasars. Following the 
method of LaMassa et al. (20151, we deredden the de¬ 
composed quasar spectrurn in the galaxy-like epoch to 
match the continuum in the quasar-like epoch to assess 
whether the broad emissi on line variability can be also be 
explained by extinction, 
that while the dimming o: 


LaMassa et al.|(2015) reported 
the quasar continuum emission 
in J0159-I-0033 (recovered in our search) can reasonably 
be modeled as due to an increase in dust extinction, the 
changes in the broad emission lines are poorly fit by the 
same extinction model. While this analysis cannot be 
performed on J0126—0839 since the quasar continuum 
dimmed completely and is not observed in the galaxy- 
like spectral epoch, we perform a similar analysis on the 
CL quasar J2336-I-0017. 

We deredden the decomposed quasar spectrum from 
the galaxy-like epoch of J2336-f0017, and fit an E{B—V) 
value for which the continuum of the dereddened spec¬ 
trum best matches the decomposed quasar spectrum 
from the earlier, quasar-like epoch of this object. This 
fitted value oi E(B — V) = 0.43 is determined by min¬ 
imizing the between the two spectra (incorporating 
all uncertainties) in the wavelength regions outside the 
H a and Bd wavel e ngth w indows discussed in Section 3.2. 
A ICardelli et al.| (|1989|) reddening law for the dust ex- 
tiriction in the host galaxy is assumed, with Ry = 3.1. 
Figure 4 compares the best-fit dereddened quasar spec¬ 
trum from the later epoch in the Ha region to that 
from the decomposed quasar spectru m from the ear¬ 
lier ep och. Similar to the findings of [LaMassa et al.| 
(20151 for J0159-I-0033, it is clear that if the dimming of 
the continuum emission between the two spectral epochs 
of J2336-I-0017 is caused purely by dust extinction, the 
change in extinction required is not consistent with the 
observed changes in the Ha emission. Furthermore, the 
profile of the broad Ha component in Figure 4 broad¬ 
ens between the two epochs of spectra. In an extinc¬ 
tion scenario, broadening of the Ha profile implies that 
the emission from the outer, lower-velocity regions of the 
broad-line region is attenuated more than emission from 
the inner portions. However, given that the quasar con¬ 
tinuum from the central accretion disk is also obscured, 
such a configuration of the obscuring material is unlikely. 
These results thus disfavor a dust extinction origin for 
the CL quasar behavior in J2336-I-0017. 


4.2. Narrow Emission Line Properties 


Merloni et al. ([2015 ) suggested that the CL quasar 
J0159-|-Ul)33 may be a I DE that was observed serendip- 
itously during and after a luminous flare, producing the 
quasar-like and galaxy-like spectra, respectively. This 
scenario is supported by their analysis of the long-term 
optical light curve of this CL quasar, which appears to 
show temporal evolution consistent with TDEs. How¬ 
ever, they also present several issues with this interpre¬ 
tation stemming from the strong broad and narrow emis¬ 
sion lines observed in the SDSS spectrum. Specifically, 
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Figure 5. BPT diagram of the three changing-look quasars in our 
sample, based on the emission line ratios measured in their latest 
SDSS spectrum. The line ratios for SDSS J015957.64-1-003310.5 
(red circle), SDSS J012648.08-083948.0 (green star), SDSS 
J233602.984-001728.7 (blue diamond), and their Ifj uncertainties 
are shown along with all emission line galaxies in SDSS-III DR12 
(blue cont ours) for comparison. T he BPT diagram classification 
schem es of|Kauffmann et al.|(|2003|l (dashed line) and|Kewley et al.| 
lf 26 oi]l (dotted iinej are snowu. i ne changing-look quasars appear 
to exhibit emission line ratios that are consistent with AGN-like or 
composite (AGN and stellar) ionizing continuum emission rather 
than powered purely by star formation alone. 


the mass of gas in the broad line region inferred by |Mer-| 


loni et al. (20151 from the spectrum (of order >100 Mq) 
is signihcantly more than could be provided from TDE 
debris. Furthermore, the possible TDE in this galaxy is 
unlikely to have ionized the gas producing the observed 
narrow lines, since the light travel tim e to the narrow 


line region (distance s of kpc scales, e.g. Liu et al(] 2013 
Hainline et al.||2013) is ^10^“'* years. Interestingly, the 


narrow line ratios of this CL quasar are not consistent 
with stellar photoionization, leaving a long-lived AGN 
as the most plausible power source. Here, we examine 
the narrow line properties for the three CL quasars in 
our sample, and compare them to those from AGN and 
TDE spectra. 

Our eigenspectra-based spectral decomposition 
method in Section 3.1 does not allow a clean separation 
of the narrow emission lines from the underlying stellar 
and quasar continuum emission. We therefore utilize the 
stellar pop ulation and emi s sion l ine fits to SDSS DR12 


spectra by Thomas et al. J 2013|^ examin e our CL 


quasars on a classic HHi' ( Baldwin et al.|1981 1 diagram. 
These fits are performed using the Gas and Abso rption 


Line Fitting (G AND ALE) software described in Sarzi 


et al. (2006), and include Gaussian fits to a variety of 


emission lines to obtain line fluxes. We specifically use 
the [N II] A6583 to Ha and [O HI] A5007 to H/3 narrow 


line ratios measured for the latest SDSS spectrum of 
our GL quasars (during the fainter, galaxy-like state), 
as shown in the BPT diagram in Figure 5. Figure 5 also 
indicates these line ratios for all galaxies in SDSS-III 
DR12 that have all four emission lines detected at >3cr 


www.sdss.org/drl2/spectro/galaxy_portsmouth/#kinematics 



Figure 6. The light curve of the continuum luminosity in the 
decomposed quasar spectra of SDSS J233602.98-1-001728.7. The 
best-fit power-law model is shown as a black solid line, and the 
best-fit power-law model with spectral index fixed to the —5/3 
value expected from tidal disruption events is shown as a black 
dashed line. 


sig nificance, as well as the BP T classification schem e 
of Kauffmann et al. (2003) and Kewley et al. (|2001 ), 
which distinguishes galaxies with emission lines ionized 
by an AGN-like, stellar-like, and composite (AGN and 
stellar) conti nuum. These l i ne ra tios for J0159-I-0033 
published by Thomas et al. (2013) are consistent with 
those inde pendently determined by [LaMassa et al.| 
(2015) and Merloni et al. (2015) for the latest SDSS 
spectrum, h igure 5 demonstrates that the emission lines 
in our sample of CL quasars are ionized at least in part 
by an AGN-like continuum, disfavoring a TDE scenario. 
For J0126—0839, the line ratios appear to be composite, 
which may be due to strong star-formation in its host 
galaxy. 

We also examine the fluxes of the narrow emission 
lines, focusing on the strongest narrow line, [ O HI] A5007. 
Lin e flux es reported for [O HI] A5007 from |Thomas~^ 


al. (2013) for the three CL quasars in our sample range 
from (1.2 - 1.7) x 10“^® erg s“^ cm“^, which corre¬ 
spond to luminosities of (1.7 - 5.3) x 10“^^ erg s“^ in 
the rest-frame. This is in strong contrast to narrow line 
emission observed in spectra of UV/optical TDEs, which 
sho w no or significantly fainter [O HI] A5007 emission 
'e.g. Gezari et al. 2006 2009[ 2012[ Holoien et al. 2014 


Chornock et al.|[20Td ), further disfavoring a 'I'DE expla¬ 


nation for GL quasars. The faint [O HI] A5007 lines 
detected in these TDEs also have ratios relative to other 
lines that are consistent with star formation rather than 
AGN photoionization; this may not be surprising given 
the long light-travel times (10^“^ years) for the ioniz¬ 
ing TDE continuum to reach the narrow line region gas. 
Thus, the combination of timescale arguments, narrow 
emission line luminosities, and emission line ratios in our 
sample of CL quasars all suggest that CL quasars are 
linked to quasar activity rather than TDEs. 


4.3. The Nuclear Light Curve of J2336+0017 

The CL quasar J2336-I-0017 in our sample has a total 
of five epochs of spectra, including four early epochs in 
its quasar-like state over approximately 2 years in the ob¬ 
served frame, and one later epoch in its galaxy-like state. 
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Although we stacked the four early epochs of spectra in 
our earlier spectral analysis since the broad emission lines 
did not show noticeable evolution, here we decompose 
separately all five epochs of spectra to study the time- 
evolution of the quasar continuum emission. A 
temporal evolution of the decaying continuum luminosity 
is often taken to be an observational signature of TDEs, 
since it is the theoretical rate at which the tidal debris 


is expected to fall back towards the SMBH 1 

Rees 

1988 

Lodato et al.|2009 

Strubbe & Quataert|2009 

1, alt 

lough 


rate (Lodato & Rossi 2011). To test the TDE explana¬ 
tion further, we compare the quasar light curve from the 
five epochs of spectra to the expected TDE luminosity 
evolution. 

Using the same decomposition method as described in 
Section 3.1, we measure the continuum luminosity at 
5100A AL 5100 of the decomposed quasar spectrum for 
the five spectroscopic epochs of J2336-I-00172. Figure 6 
presents the light curve, along with the best-fit power- 
law of i-i i3±o.20^ best-fit TDE model. 

The uncertainty on this power-law index is calculated by 
10^ Monte Carlo resamplings of each point on the light 
curve from a Gaussian based on their Icr uncertainties, 
to produce 10^ resampled light curves. Each resampled 
light curve is refit, and the quoted Icr uncertainty on 
the power-law index is the Icr spread from the resampled 
light curves. The light curve’s best-fit power-law index of 
— 1.13 ±0.20 is statically shallower than the —5/3 model 
predicted for TDEs, although this only weakly disfavors 
the TDE model due to the sparse light curve sampling 
and uncertainties in the models for temporal evolution 
of TDE luminosity. 


4.4. A Type Iln Supernovae Origin? 

Type Iln supernovae (SNe) frequently display broad 
H Balmer emission and blue UV continuum emission in 
their optical/UV spectra, which can appear qualitatively 
similar to AGN spectra du ring certain phases of their 
evolution (|Filippenko 1989|). Thus, it may be possible 
that our CL quasars are actually serendipitously discov¬ 
eries of Type Iln SNe. This possibility is ruled out for 
two of our CL quasars (J0159±0033 and J2336±0017) 
because broad Ha emission is still detected in their later 
galaxy-like spectra approximately 7 years later in the 
rest-frame. Furthermore, the narrow line ratios for these 
two objects are clearly AGN-like (see Figure 5), indicat¬ 
ing the presence of nuclear activity. For J0126—0839, no 
broad emission lines or quasar continuum emission is de¬ 
tected in the later galaxy-like spectral epoch, leaving this 
dimming timescale poorly constrained. Furthermore, the 
narrow line ratios for this object fall in the ‘composite’ 
rather than the AGN region of the relevant diagnostic 
line ratio diagrams, making a Type Iln SN explanation 
for this object somewhat more plausible at first glance. 
However, the Ha luminosity of J0126-0839 in the bright 
state (10^^ ° erg s“^ from Table 1) is approximately 1-2 
orders of magnitude more than the peak Ha luminosi¬ 
ties normally observed in Type H n supernovae (e.g. see 
Figure 14 of Taddia et al.||2013 ). The strong Ha emis¬ 
sion observed in Type Hn supernovae emerges several 
months after their light curve peaks, when the contin¬ 
uum emission has dimmed by several magnitudes and is 


red ( Filippenko||1997 ). In contrast, the strong Ha emis¬ 
sion observed in the bright-state spectrum of J0126-0839 
is coincident with highly-luminous and blue continuum 
emission. P Cygni profiles and other asymmetries of¬ 
ten observed in the Ha line of Type Hn SNe are also 
not seen in our GL quasars. For these reasons, we dis¬ 
favor the Type Hn supernovae scenario for J0126—0839, 
but this possibility cannot be conclusively ruled out be¬ 
cause the nature of our search for GL quasars using a 
large data set may preferentially uncover such rare phe¬ 
nomena. Follow-up X-ray and radio observations of the 
nucleus of J0126—0839 may be useful to confirm its AGN 
nature. 


4.5. Infall Timescales 

Our investigations above generally disfavor TDEs, 
Type Hn SNe, and dust obscuration as the origin of GL 
quasars, and instead suggest that this phenomenon is due 
to intrinsic dimming of the quasar emission from rapidly 
decreasing accretion rates. Here, we assess whether the 
±10 year transition timescales we observe for our sample 
of GL quasars are consistent with the infall timescale of 
gas in the radiatio n-pressure dominated inne r regions of 
Shakura-Sunyaev (Shakura & Sunyaev||1973 ) thin accre¬ 
tion disks; this is also the timescale on which changes in 
the accretion rate are reflected in change s in the contin¬ 
uum l uminosity. Using Equation 5 from |LaMassa et al.| 
(20151, we find that the infall timescale for our sample of 
changing look quasars are approximately 42, 38, and 868 
years for J0159±0033, J0126-0839, and J2336±0017, re¬ 
spectively. Although these infall timescales are longer 
than the transition timescales we observed, we note that 
the transition timescales observed in our CL quasars are 
lower limits since the multi-epoch spectra do not encom¬ 
pass the full transition (e.g., it is likely that the CL 
quasars have begun dimming before the fir st spectral 
epoch). Furthermore, LaMassa et al. (2015) also point 
out that magneto-hydrodynamic simulations of quasar 
accretion flows have suggested that the infall timescale 
may be a factor of a few sho rter than these analytical 
estimates (Krolik et al. 2005), and much closer to the 
observed Ai(j year timescales. However, the 868 year in¬ 
fall timescale estimated for J2336±0017 is problematic 
in this interpretation, and may indicate that other pro¬ 
cesses such as thermal or dynamical instabilities may be 
present in the accretion flows of CL quasars. 


5. CONCLUSIONS 

The discovery of CL quasars presents a new oppor¬ 
tunity to study the nuclear environment and structure 
of quasars, once the origin of this phenomenon is un¬ 
derstood. To provide a larger sample of these objects, 
we performed an archival spectroscopic search in SDSS, 
yielding three CL quasars, including two new cases. Us¬ 
ing this sample, we investigate the detailed properties of 
their quasar continuum emission and broad and narrow 
emission lines, with the goal of attempting to discrim¬ 
inate between various possibilities for the origin of this 
phenomenon. The primary results of our investigation 
can be summarized as below: 


• The three CL quasars in our sample appear to 
show similar properties: they are luminous (Lboi ~ 
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1044 o- 44.5 gj.g g-i^ quasars at relatively low red- 
shifts (z ~ 0.2 — 0.3) that display strong dimming 
of the quasar continuum and the broad H Balmer 
emission lines over timescales of approximately 5 
to 7 years in the rest-frame. 

• In the two CL quasars for which the broad emis¬ 
sion lines are detectable in both spectral epochs, 
the decrease in broad line luminosity coincides 
with broadening of the broad line widths, such 
that the derived black hole masses are preserved. 
This is consistent with a rapid change in their Ed¬ 
dington ratios, which decrease from Lboi/LEdd 
0.03 — 0.005 until the broad H Balmer lines have 
dimmed significantly or disappeared. In one CL 
quasar (J0126-0839), all quasar continuum emis¬ 
sion and broad emission lines have disappeared be¬ 
low the SDSS detection limit. 


• Changes in dust extinction required to match the 
dimming in the quasar continuum cannot account 
for the changes in the broad emission lines in ei¬ 
ther of our two CL quasars for which this anal¬ 
ysis was possible, disfavoring an extinction origin 
for this phenomenon. Narrow emission line diag¬ 
nostics show that our CL quasars all have lumi¬ 
nous narrow lines with line ratios consistent with 
at least partially AGN-like ionizing emission. We 
argue that these narrow line properties favor a sce¬ 
nario in which the quasar continuum dims intrinsi¬ 
cally over a TDE origin for this phenomenon. 


If the intrinsic dimming of the quasar emission fa¬ 
vored by our analysis is due to draining of the underlying 
quasar accretion disk, the CL quasar phenomenon will 
provide a unique new laboratory to study the accretion 
flow and nuclear environment in luminous AGNs. Long¬ 
term spectroscopic and multi-wavelength monitoring of 
the currently-known changing-look quasars can help fur¬ 
ther elucidate the origin of CL quasar transitions. For 
example, observations of the return of any CL quasar to 
a bright quasar-like state would provide additional con¬ 
straints on the physical mechanism of this phenomenon 
as well as estimates of its duty cycle. Current multi¬ 
object spectroscopic programs and time-domain imaging 
surveys are well-poised to serendipitously discover many 
more CL quasars. For example, the Time-D omain Spec¬ 
troscopic Survey (Morganson et al. 2015) in SDSS-IV 
will provide repeat spectroscopy ot several thousand low- 
redshift quasars, while the Pan-STARRS 37r survey has 
repeated imaged 30,000 deg^ of sky (including the SDSS 
imaging footprint). With future instruments and sur- 


(DESI, Levi et al.J 

2013) and the Large Svnoptic Survev 

Telescope (LSS'L [ 

Lvezic et al. 

2008 

1, discovery of such 
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Figure 7. SDSS and APO 3.5m rest-frame spectra of SDSS J132457.29-1-480241.2. The straight red line in the SDSS spectrum of this 
object at MJD = 56805 is directly from the SDSS spectrum file, and indicates that this spectrum was not obtained or reduced correctly. 
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APPENDIX 

A. THE CL QUASAR CANDIDATE SDSS J132457.29-^480241.2 

In addition to the three CL quasars yielded by our systematic search, we identified one additional possible CL quasar 
in our visual inspection (SDSS J132457.29-I-48024I.2) that is classified as CLASS = ‘QSD’ in its earlier SDSS spectrum, 
and CLASS = 'GALAXY’ in its latest SDSS spectrum (shown in Eigure 7). However, the red side of the latest SDSS 
spectrum is corrupted beyond SOOOA, although the blue side appears to show the disappearance of broad H/3 and Hy 
emission, along with dimming of the quasar continuum, similar to CL quasars. The fiber (527) of the latest SDSS 
spectrum is known to be affected by columns of bad pixels in the red camera of the BOSS spectrograph, which is likely 
to be the cause of the corrupted spectrum (and affects spectra from neighboring fibers on this and other plates). 
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Nevertheless, to verify whether this object is indeed a CL quasar, we obtained additional optical long-slit spectra 
using the Dual Imaging Spectrograph on the Astrophysical Research Consortium 3.5m telescope at Apache Point 
Observatory, with wavelength coverage of A ~ 3400 — 9200 A. Three 15 minute exposures where taken on 14 January 
2015 UT, at a spectral resolution of R~800 using the B400/R300 grating settings, and a 1.5" slit. The seeing was 1.8" 
on this night, and the observations were obtained at airmass of approximately 1.04. Spectra of the spectrophotometric 
standard star Feige 34 were also obtained for flux-calibration and removal of atmospheric absorption, and HeNeAr 
lamps were used to obtain a wavelength solution. These spectra were bias and flat-held corrected, wavelength- and 
flux-calibrated, and corrected for atmospheric extinction using standard IRAF procedures. 

The calibrated APO 3.5m spectrum is displayed in the lower panel of Figure 7, which shows that although the quasar 
continuum, broad Ha, and broad H/3 emission has dimmed, they remain prominent, and thus we do not include this 
object in our CL quasar sample. Although the disagreement between the APO 3.5m spectrum and the latest SDSS 
spectrum is likely due to the known column of bad pixels, a scenario in which this object is a CL quasar that dimmed 
in the latest SDSS spectrum to a galaxy-like state, then rebrightened back to a quasar-like state in the APO 3.5m 
spectrum cannot be ruled out. Given the short period between the latest SDSS spectrum and APO spectrum (182 
days in the rest-frame), and the known issues with the latest SDSS spectrum in the red side, rebrightening of a CL 
quasar appears unlikely, but additional spectral monitoring may be warranted. 


B. REJECTED CL QUASAR CANDIDATES 

Our systematic search for CL quasars based on multi-epoch SDSS spectra (outlined in Section 2.1) yielded 117 
candidates for which the SDSS spectral pipeline classification switched between CLASS = ‘GALAXY’ and CLASS = 
‘QSO’ (or vice versa) in repeat spectra. Visual inspection of these 117 candidates resulted in the 3 confident CL 
quasar candidates, while the remaining 114 were rejected because they did not display the disappearance of any broad 
line and continuum emission characteristic of CL quasars. Figure 8 shows example repeat spectra of some of these 
candidates that were rejected in our visual inspection, which includes 24 objects in which broad Ha emission was 
outside the spectral range in one spectral epoch but not the other. Figure 8 also shows the repeat spectra of one object 
(SDSS J032332.83—000740.3) for which the fiber was positionally offset between the two spectral epochs. A positional 
offset of the fiber (e.g. off the nuclear region) can artificially cause changing-look behavior, although this would also 
cause significant artificial variability in the host-galaxy emission. 


C. TRANSITION TYPES 


Our search criteria for CL quasars should allow us to And cases in which objects spectroscopically transformed from 
quasar-like to galaxy-like, as well as the reverse. However, in our visual inspection of the CL quasar candidates, we And 
no confident cases of galaxy-like to quasar-like CL qua sars. Although this reverse transition (galaxy-like to quasar-like) 
has been observed previously in Seyfert galaxies (e.g., Goodrich|1995 Shappee et al.|2014 ), these Seyferts are at lower 
redshifts and lower luminosities than our CL quasarsl It is thus unclear whether the lack of reverse transitions in 
our sample is due to selection effects, a random result stemming from our small sample size, or a true feature of CL 
quasars. 

We first consider potential selection effects that may cause us to select CL quasars of only the quasar-like to galaxy¬ 
like variety. For example, if SDSS systematically targeted quasars for repeat spectroscopy in greater numbers than 
galaxies, then naturally we would find more quasar-like to galaxy-like transitions. To test this possibility, we identify 
the number of repeat spectra in SDSS DR12 of objects with CLASS = ‘GALAXY’ and objects with CLASS = ‘QSO’, 
using the same search criteria described in Section 2.1. There are 11,438 repeat spectra of 8,865 unique galaxies, and 
7,109 repeat spectra of 5,990 quasars (some objects had more than two epochs of spectra). This test shows that our 
sample of only quasar-like to galaxy-like CL quasars is not simply due to a lack of repeat galaxy spectra in SDSS. 
However, the myriad other potential selection effects, such as those stemming from the disparate SDSS galaxy and 
quasar completeness magnitude limits, the quasar and galaxy targeting methods in SDSS, and differences in the galaxy 
populations between those that host quasars and those targeted by SDSS spectroscopy make it difficult to draw robust 
conclusions on selection effects without more careful considerations. 

We next consider whether our yield of three CL quasars of the quasar-like to galaxy-like variety is simply a random 
result due to the small number of CL quasars thus far discovered. Under the simple assumption that transitions of 
either variety are equally likely both intrinsically and observationally (i.e. neglecting selection effects), the probability 
of observing three quasar-like to galaxy-like transition and no cases of the reverse is approximately 0.5^ = 12.5% from 
the binomial distribution. In the absence of selection effects, it thus appears unlikely for our search to randomly yield 
only quasar-like to galaxy-like CL quasars if the reverse transition is equally likely, though the statistics are not yet 
compelling. Recent searches for CL quasars using alternative approach es based on photomet ric light curves have now 
uncovered objects interpreted as galaxy-like to quasar-like CL quasars (MacLeod et al.|[2015). 
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Figure 8. Example repeat SDSS spectra of the CL quasar candidates that were rejected in our visual inspection. The MJD and SDSS 
pipeline classification of each spectral epoch are shown. These rejected candidates include objects in which the broad Ha emission is 
outside the spectral range in one epoch but not the other (e.g., SDSS J230817.30—002128.8 and SDSS J023359.92+004012.7). For one 
object (SDSS J032332.83—000740.3), the fiber was positionally offset, causing the observed artificial variability in the host-galaxy emission. 
Aside from this object and the three CL quasars, all other candidates were rejected because our visual inspection of their repeat spectra 
did not reveal dramatic disappearance or appearance of broad emission lines. 







